Holder Mean applied to Anderson localization 
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The phase diagram of correlated, disordered electron systems is calculated within dynamical 
mean-field theory using the Holder mean local density of states. A critical disorder strength is 
determined in the Anderson-Falicov-Kimball model and the arithmetically and the geometrically 
averages are found to be just particular means used respectively to detect or not the Anderson 
localization. Correlated metal, Mott insulator and Anderson insulator phases, as well as coexistence 
and crossover regimes are analyzed in this new perspective. 
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I. INTRODUCTION 

The properties of materials are strongly influenced 
by electronic interaction and randomness. In particu- 
lar, Coulomb correlations and disorder are both driv- 
ing forces behind metal-insulator transitions (MIT) con- 
nected with the localization and derealization of par- 
ticles. The Mott-Hubbard MIT is caused by Coulomb 
correlations (electronic repulsion) in the pure system [l[ . 
The Anderson MIT, also referred to as Anderson local- 
ization, was established in [2| and is the basis of the 
theory of localization of electrons in disordered systems. 
It is due to coherent backscattcring from randomly dis- 
tributed impurities in a system without interaction Q. 
It is therefore a challenge to investigate quantum models 
where both aspects are simultaneously present. 

The Mott-Hubbard MIT is characterized by opening a 
gap in the density of states at the Fermi level. At the 
Anderson localization, the character of the spectrum at 
the Fermi level changes from a continuous one to a dense 
discrete one. It is plausible that both MITs could be de- 
tected by knowing the local density of states (LDOS), as 
it discriminates between a metal and an insulator, which 
is driven by correlations and disorder Q • In many recent 
works [a, |y, lZ|, the change of the whole LDOS distribu- 
tion function at the Anderson transition were explored 
for the Bethe and for the simple cubic lattices, where 
many other effects, like interaction, percolation and bi- 
nary alloy disorder, were included. 

The theoretical investigation of disordered systems 
requires the use of probability distribution functions 
(PDFs) for the random quantities of interest and one 
is usually interested in typical values of these quantities 
which are mathematically given by the most probable 
value of the PDF. In many cases the complete PDF is not 
known, i.e., only limited information about the system 
provided by certain averages (moments or cumulants) is 
available. In this situation it is of great importance to 
choose the most informative average of a random vari- 
able. For examples, when a disordered system is near 
the Anderson MIT, most of the electronic quantities fluc- 



tuate strongly and the corresponding PDFs possess long 
tails Q and at the Anderson MIT the corresponding mo- 
ments might not even exist. This is well illustrated by 
the LDOS of the system. The arithmetic mean of this 
random one-particle quantity does not resemble its typ- 
ical value at all. In particular, it is noncritical at the 
Anderson transition and hence cannot help to detect the 
localization transition. By contrast, it was shown 0, Q 
that the geometric mean gives a better approximation of 
the averaged value of the LDOS, as it vanishes at a criti- 
cal strength of the disorder and hence provides an explicit 
criterion for Anderson localization [1, Gil [III> [j3 • Besides 
that, the work [l3| suggests that by considering the An- 
derson localization, the moments of the distribution of 
the resolvent kernel are of essential importance. 

Theoretical descriptions of the MIT have to be non- 
perturbative if no long-range order exists on either side 
of the transition. A nonperturbative framework to inves- 
tigate the Mott-Hubbard MIT in lattice electrons with a 
local interaction and disorder is given by the dynamical 
mean-field theory (DMFT) pj]. 

In this paper, we investigate the Anderson-Falicov- 
Kimball model. The pure Falicov-Kimball model de- 
scribes two species of particles, mobile and immobile, 
which interact with each other when both are on the 
same lattice site. The Falicov-Kimball model captures 
some aspects of the Mott-Hubbard MIT, i.e., upon in- 
creasing the interaction the LDOS for mobile particles 
splits into two subbands opening a correlation gap at the 
Fermi level. In the Anderson-Falicov-Kimball model the 
mobile particles are disturbed by a local random poten- 
tial. The model is solved within the DMFT framework. 

Anyway, there is no reason to believe that just the 
geometric mean can offer a good approximation for the 
averaged LDOS. Actually, if we consider the generalized 
mean we can show that the averaged LDOS can vanish in 
the band center at a critical strength of the disorder for a 
wide variety of averages. Here, we use the Holder mean in 
order to analyze how the averaged LDOS depends on each 
Holder parameter that is used and whether, in particu- 
lar, the geometric mean can offer a better approximation 



than the arithmetic one. 

A generalized mean, also known as Holder mean, is an 
abstraction of the Pythagorean means including arith- 
metic, geometric and harmonic means. If q is a non-zero 
real number, we can define the generalized mean with 
exponent q of the positive real numbers x±, . . . , x n as 
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As special cases, we can find, for example, 
the minimum (limq^-oo M q (x)), the geometric mean 
(liniq^o M q (x)), the arithmetic mean (M%(x)) and the 
maximum (lim g _»oo M q (x)). A power mean emphasizes 
small (big) values for small (big) q. 

In section [II] we present the Anderson- Falikov-Kimball 
model and its solution using the DMFT [9J. We also 
present the relation between the disorder strength and 
the Coulombian repulsion for the linearized DMFT with 
g-averaging. In section IIIII we present the numerical 
results concerning the ground-state phase diagram ob- 
tained using the Holder mean and compare them with 
the results of previous works 4. 9. 11]. Finally in section 
IIVI we present our conclusions and final remarks. 

II. ANDERSON- FALICOV-KIMB ALL MODEL 



one electron occupy the same site. The average num- 
ber of electrons (ions) on site i is denoted as n e — cfa 
(rif = / 4 + /i).The energy e$ is a random, independent 
variable, describing the local disorder disturbing the mo- 
tion of electrons. For simplicity, it was assumed that just 
mobile particles are subjected to the structural disorder 

The number of ions and electrons is considered inde- 
pendent of each other and fixed. If there is no long- 
range order, the position of the ions on a lattice is ran- 
dom, what introduces additional disorder apart from that 
given by the term in the Hamiltonian We consider 
that the occupation rif on the ith site has probability p 
(0<p<l). 

In the U term in the Hamiltonian ^ one has to take 
the quantum-mechanical average over a given quantum 
state of the / particles, what does not transform the ex- 
tended states into the localized ones. In the term one 
has to average the quantum-mechanical expectation val- 
ues over different realizations of e*, what can lead to the 
Anderson localization Q. 

In the pure Falicov-Kimball model, if n e + n/ = 1, 
the Fermi energy for electrons is inside of the correlation 
(Mott) gap opened by increasing the interaction [9( . How 
the disorder changes this gap and how the g-averaging 
influences the presence of localized states are the subjects 
of the present study. 



The Model 



B. Dynamical mean-field theory 



The Falicov-Kimball model is, presently, the simplest 
model to study metal-insulator transitions in mixed va- 
lence compounds of rare earth and transition metal ox- 
ides, ordering in mixed valence systems, order-disorder 
transitions in binary alloys, itinerant magnetism, and 
crystallization. Recently, it was also applied to study the 
possibility of electronic ferroelectricity in mixed-valence 
compounds, and also of the phase diagram of metal 
ammonia solutions [l4| . In its most simplified version, 
namely the static model, it consists in assuming that in 
the system exist two species of spinless fermions: one of 
them possess infinite mass and hence does not move while 
the other one is free to move [To| . 

In the Anderson- Falicov-Kimball model the mobile 
particles are disturbed by a local random potential, giv- 
ing rise to a competition between interaction and disorder 
yielding stabilization of metalicity. This model is defined 
by the following Hamiltonian: 



cfci 



(2) 



where (q) and ff' (ft) are, respectively, the cre- 
ation (annihilation) operators for the mobile and immo- 
bile fermions (electrons and ions, respectively) at a lat- 
tice site i, t is the electron transfer integral connecting 
states localized on nearest neighbor sites and U is the 
Coulombian repulsion that operates when one ion and 



The Anderson- Falicov-Kimball model, where the inter- 
action and disorder are local, is solved within the DMFT. 
The formalism that is used is the same from Q , where ui 
denotes the energy, the number of the immobile particles 
is conserved being zero or one and the chemical potential 
fi is introduced only for the mobile subsystem. 

Using the hybridization function r){uS), which is a dy- 
namical mean field (molecular field) describing the cou- 
pling of a selected lattice site with the rest of the system 



pimg ot a selected lattice site witn 
[III], the ei-dependent LDOS is @ 
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where G(u>, eA is the local e^-dependent Green function 
described in [9(. 

From the e^-dependent LDOS, we introduce the q- 
averaged LDOS 
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where the subscript q stands for the chosen average. 

The lattice, translationally invariant, Green function 
is given by the corresponding Hilbert transform 



G(u) = / dio 



(5) 
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The self-consistent DMFT equations are closed 
through the Hilbert transform G(w) = J deN (e)/{—e + 
7](lu) + 1/G(w)], where N (e) is the non-interacting den- 
sity of states. Using the semi-elliptic density of states for 
the Bethe lattice, that is N (e) = 4^1 - A{t/W) 2 / {ttW) 
fTst ] , the Anderson- Falicov-Kimball model can be exactly 
solved. In this case, we can find ri(ui) = W 2 G{us)/ld> Q. 
We considered ti an independent random variables char- 
acterized by a probability function P(ti). We assume 
here a box model, i.e., P(ej) = <&(A/2 — e^/A, with $ as 
the step function. The parameter A is a measure for the 
disorder strength. If we consider the real and imaginary 
parts of rjiuj) we can use r(uj) +is(uj) — G ^ , what leads 
to 



a? + s 2 + {U/2f 



it [a? + s 2 + (17/2) 



212 



U 2 a 2 



(6) 



where on = lu — e$ — r. 

The chemical potential p = U/2, corresponding to a 
half- filled band (i.e., n e = 1/2), andp =1/2 are assumed 
in this paper. W = 1 sets the energy units. 



C. Linearized Dynamical Mean-Field Theory 

At the MIT, the LDOS vanishes in the band center 
and is arbitrarily small in the vicinity of the MIT (but 
on the metallic side). The ground-state properties in 
the half-filled band case are solely determined by the 
quantum states in the band center (u> = 0). In this 
sense, by linearizing the DMFT equations we can de- 
termine the transition points on the phase diagram. Due 
to the symmetry of Pq{oS), in the band center we have 
that G(0) = —iirp q {ijj) and is purely imaginary [9(, what 
leads to the recursive relation 7/ ra+1 )(0) = — i7rp"(0)/16. 
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FIG. 1: Ground-state phase diagram for electrons in a band 
center determined by using the g-average. Dots are deter- 
mined from the numerical solution of the DMFT equations. 
Solid lines are obtained by numerical integration from the 
linearized DMFT. 



The left hand side in the (n + l)th iteration step is given 
by the result from the (n)th iteration step. Using Eq. 
([6|) and expanding it with respect to small p™ (0), we find 
that the recursive relations within the linearized DMFT 
with g-averaging are 
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where 
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At the boundary curves between metallic and insu- 
lating solutions the recursions are constant, p™ +1 (0) = 
Pq(0), because in a metallic phase the recursions are in- 
creasing, whereas in the insulating phase they are de- 
creasing. This observation leads directly to the ex- 
act (within DMFT) equations determining the curves 
A q = A q (U), i.e., 
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for the linearized DMFT with ^-averaging. 



III. RESULTS 
A. The critical disorder strength 

In the limit U = 0, the equation (J9j) can be analytically 
solved and we can find the critical A as a function of q 
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FIG. 2: Ag as a function of q for (7 = 0. The point represents 
the exact result for the cubic lattice corresponding to q = 



0.011424 and 



1.375. 
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We plot this function in Fig. [2l As A^(0) is real, finite 
and positive, we must have q < 1/2. 

This shows that the Anderson localizations can be de- 
tected when we consider a generalized mean with q < 
1/2. In this context, one can easily note that the geo- 
metrical and the arithmetical averages (q = and q = 1, 
respectively) are just particular means, where, respec- 
tively, the Anderson localization can either be detected 
or not. 

The numerical integration of equation (|9|) is shown as 
solid curves in Fig. Q]for different q's. We can clearly see 
that for values of q greater than 0.5, A always increases 
with U, what characterizes the Mott-Hubbard MIT. For 
values of q smaller than 0.5 we can always identify three 
regions with respect to the states at the Fermi level: ex- 
tended gapless phase, localized gapless phase and gap 
phase [3]. 

As already pointed out in 9] , depending on the spec- 
tral properties of the Anderson-Falicov-Kimball model, 
we distinguish three different regimes: weak, intermedi- 
ate and strong interaction regime. In the weak interac- 
tion regime the Mott gap is not open; in the intermediate 
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FIG. 3: Averaged local density of states at q = 0.5 for U — 0.3 
(upper panel) and U = 0.9 (lower panel) for different disorder 
strengths A (0, 0.4, 0.8, 1.2, 1.6, 2.0 and 2.4). 
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FIG. 4: Averaged local density of states at U — 0.3, for A = 
1.4 and q = -0.25, 0, 0.25, 0.5, 0.75, 1.0 and 2.0. As q 
decreases A q (0) is smaller. 



interaction regime, the Mott gap is open at A = and in 
the strong interaction regime the Mott gap, determined 
within the DMFT framework with q less than 0.5, is al- 
ways open even at large disorder. 

The Mott and Anderson insulators are continuously 
connected. Hence, by changing U and A it is possible 
to move from one type of insulator to the other without 
crossing the metallic phase. This occurs because the An- 
derson MIT(E7 = 0) and the Mott-Hubbard MIT (A = 0) 
are not associated with a symmetry breaking [9j . 

When we look at the curves obtained for q smaller than 
0.5, where the Anderson localization can be detected, we 
see (Fig. [1]) that the curves obtained for greater values 
of q enclose those ones for smaller values of q. If we look 
at the curve corresponding to the smallest value of q we 
see that in order to obtain the Anderson localization, one 
needs to consider a small range of values of U and A; in 
the other extreme, for values of q close to 0.5 the range 
of values of U that can be considered is large, but one 
needs extremely high values of A to detect it. 



B. The LDOS 

In order to understand how the numerical solutions of 
the curves of Fig. [T] are obtained one needs to observe at 
first the behavior of the average local density of states for 
different values of U as the disorder strength A is varied. 
This is made considering equations ([3])- (J5J), from which 
we can obtain the stable value of pq(u>) as a function of 
u>, as can be seen in Fig. [3l 

In all simulations we considered that the initial value of 
p(tjj, £j) is a uniform distribution with bandwidth W = 2t 
and then we determined G(u>) in order to obtain rj(u>) and 
finally the new values of p(ui, ej). This loop is performed 
until we find the stable configuration for p q {oj). 

If we consider fixed values of U and A we can observe 
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FIG. 5: Averaged local density of states at U = 0.9 for q = 
-0.25, 0, 0.25, 0.5, 0.75, 1.0 e 2.0 and for (a) A = 0.6, (b) 0.8 
and (c) 1.2. As q decreases A q (0) is smaller. 



the dependence of the averaged LDOS on the value of 
q that is used. Some of these results are presented in 
Fig. 0] and Fig. [5] In the first one, the simulations are 
performed in a weak interaction regime (U = 0.3). In 
the second one, we consider an intermediate interaction 
regime (U = 0.9). Although the behavior of the LDOS is 
not the same in both cases, we can see that the smaller 
the value of q that we consider, the smaller is the value of 



p q (oj = 0). We can also see that the shape of the curves 
change when we compare the curves for q smaller than 
0.5 and greater than 0.5. 

A signature of the Anderson localization is the vanish- 
ing of p q (uj = 0) as we increase A. The way to obtain the 
results in Fig. [T] consists in considering a range of values 
of U for a fixed q and varying the disorder strength A 
for each value of U, in order to determine the values of 
A when p q (uj = 0) = 0. 

In the simulations, we used values of A varying of 0.2 
and determined the stable values of p q (io = 0) for each 
of them. As we use an iterative process, the values of 
p q (u> = 0) that we obtain converge always to the stable 
one. In this sense, for some values of A the stable value 
of pq[io = 0) is exactly null. In Fig. [6] we show the de- 
pendence of p q {lo = 0) on A for two different values of 
U (0.3 and 0.9) that, again refers to the weak and in- 
termediate interaction regime, respectively. For the first 
one, increasing A at fixed U and q leads to further de- 
creasing of the LDOS; just for values of q smaller than 
0.5, p q {uj — 0) can vanish for a finite value of A; for q 
greater than 0.5 this would just happen at infinity. For 
q smaller than 0.5 we found that the critical exponent 
of the curves p q (oj — 0) « (A — Ac)' 3 is always (3 = 1 




FIG. 6: Average density of states in a band center (ui = 0) 
as a function of disorder A with U = 0.3 (upper panel) and 
U = 0.9 (lower panel). The results are obtained for q = —0.25, 
0.25, 0.5, 0.75, 1.0 e 2.0. 
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which guarantees that they all belong to the same 
universality class. For the second plot, p q (u> — 0) is zero 
until a certain value is reached, which corresponds to the 
transition from the gap phase to the extended gapless 
phase. Here again just for values of q smaller than 0.5, 
p q (u> — 0) turns to vanish for a finite value of A. Per- 
forming the same process for other values of U we obtain 
easily all the simulation points of Fig. [1] for a variety of 
different values of q. 

IV. CONCLUSIONS 

In the present paper, we generalized the solutions 
of the Anderson-Falicov-Kimball model using a Holder 
mean for the averaged local density of states. In this 
context, where each mean is characterized by the Holder 
parameter q, we found that the averaged LDOS can be 
calculated for arbitrary parameters in order to detect the 
Anderson-localization. As we vary the parameter in the 
Holder mean, we emphasize different signal- values of the 
LDOS. 

Besides that, there is a critical value for this parameter, 



that is q = 1/2. For values of q greater and smaller than 
this, the dependence of the critical strength of disorder 
A on the Coulombian repulsion U is, respectively, similar 
to the one found in Q using the arithmetic mean (q = 1) 
and the geometric mean (q — 0). We showed that not 
only the geometric mean can be used to detect the An- 
derson mean, but also an infinity of other averages. One 
example occurs for the cubic lattice [ll], [l7| , where the 
exact solution A^ = 1.375 corresponds to q = 0.011424 
in our model (Fig. [2]) , which is not exactly the geometric 
mean. 

In this way, using the Holder mean, we were able to 
understand why in many previous works the geometric 
mean has offered a better approximation for the averaged 
LDOS than the arithmetic one. 
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